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Abstract 

In pQ, the authors consider the generalization Q y K of the Erdos-Renyi random graph process 
Qi, where instead of adding new edges uniformly, gives a weight of size 1 to missing edges 
between pairs of isolated vertices, and a weight of size K S [0, oo) otherwise. This can correspond 
to the linking of settlements or the spreading of an epidemic. The authors investigate tg(K), 
the critical time for the appearance of a giant component as a function of K, and prove that 
tg=(l + o(l)) J^ K , using a proper timescale. 

In this work, we show that a natural variation of the model has interesting properties. 
Define the process Q^, where a weight of size K is assigned to edges between pairs of non- 
isolated vertices, and a weight of size 1 otherwise. We prove that the asymptotical behavior 
of the giant component threshold is essentially the same for and namely tg/t^ tends to 
tt(24+^) ps 1.47 as K — > oo. However, the corresponding thresholds for connectivity satisfy 
*c Ac = max{i,.ftT} for every K > 0. Following the methods of pQ, tg is characterized as the 
singularity point to a system of differential equations, and computer simulations of both models 
agree with the analytical results as well as with the asymptotic analysis. In the process, we 
answer the following question: when does a giant component emerge in a graph process where 
edges are chosen uniformly out of all edges incident to isolated vertices, while such exist, and 
otherwise uniformly? This corresponds to the value of tg(Q), which we show to be | + 3c2 _ 1 . 

1 Introduction 



1.1 The Erdos-Renyi graph process and biased processes 

The random graph process, Gi(n), is a sequence of (2) + 1 graphs on n vertices, Q\°, . . . , gS 2 \ where 
Qi° is the edgeless graph on n vertices, and Q\ T is obtained by adding an edge to G\ T ~ l , chosen 
uniformly over all missing edges. This model was introduced by Erdos and Renyi in where it is 
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shown that for every constant C < |, the largest component of Q\ is typically of size O(logn), 
and yet for every constant C > | there is typically a single component of linear size in Gi Cn and all 
other components are of size O(logra). This single component and its evolution are referred to as 
the "giant component", and "the double jump phenomenon" respectively. Using a timescale of n/2 
edges, the threshold for the appearance of the giant component is thus t g = 1. Another classical 
result of Erdos and Renyi determines that g 1 Cnlo s n \ s typically connected for C > ^, and typically 
disconnected for C < ^. Using a timescale of ^ logn, the connectivity threshold is thus t c = 1. For 
further information on the evolution of the giant component in the random graph process, as well 
as on the threshold for connectivity, see, e.g., 

There has been extensive study on the thresholds for the appearance of a giant component 
and for connectivity in different variations of the random graph process. For instance, in a model 
suggested by Achlioptas, two random edges are chosen uniformly out of the missing edges at each 
step, out of which some algorithm A selects one to be added to the graph. For results on upper 
and lower bounds on the emerging of the giant components for various algorithms in this model, 
see @|, 0,|5|, HDI- 

In [Jj, the authors consider a random graph process on multi-graphs, where at each step an 
edge is added between a random vertex of minimal degree and a random uniformly chosen vertex. 
The authors analyze the number of vertices of degrees 0, 1, 2 along the process using the differential 
equation method for graph processes of Wormald 12 J, and show that the mentioned graph process 
becomes connected typically when the minimal degree becomes 3. 

The following generalization of the original graph process Qi, which we denote by Q^{n), was 
studied in pQ: at each step, Q y K gives a weight of size 1 to missing edges between pairs of isolated 
vertices, and a weight of size K € [0, oo) otherwise (when no isolated vertices are left, the distribu- 
tion on the missing edges becomes the uniform distribution). This can correspond to the linking 
process of n initially isolated settlements, or the spreading of an epidemic, where the probability of 
a new link is affected by whether or not one of its endpoints already has other links. The threshold 
for the appearance of a giant component in Q^, t v , becomes a continuous function of the parameter 
K, and the authors of [Q use the differential equation method to express t y as a singularity point 
to a system of coupled non-linear ordinary differential equations (ODEs). By applying methods of 
asymptotic analysis of ODEs, it is proved that t y g {K) = (1 + o(l)) -7=^, where the o(l)-term tends 
to as K — > 00. 

In this work, we show that a natural variation on the process Q y K has interesting properties. 
Consider the process Q^-, where instead of placing the weight K when one of the endpoints is 
non-isolated (as in Gk), it is placed when both endpoints are non-isolated. In other words, Q y K gives 
a weight of size K S [0, 00) to missing edges between pairs of non-isolated vertices, and a weight of 
size 1 otherwise. Notice that for K = 1, both processes are equivalent to the original graph process. 
Furthermore, for any K, both processes Q y K and apply the rule of the original graph process 
Qi once no isolated vertices are left, hence it is interesting to compare the two up to roughly that 
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point. 

Applying the methods of on Qg, we express its threshold for the appearance of a giant 
component, tg(K), as a singularity point to a system of coupled non-linear ODEs. For the special 
case K = 0, we obtain a setting where edges are added uniformly at random out of all edges incident 
to an isolated vertex, until no such vertex is left (from that point on, edges are added uniformly at 
random). In this special case we get tg (0) = | + 3e -f_ 1 ■ 

Using classical methods of asymptotic analysis of ODEs (see, e.g., [2]) we readily calculate the 
asymptotic behavior of t£(K) for K » 1, and obtain that t£(K) = (1 + o(l)) + ~Jk~ 

It follows that t^/t y g « 1.47 for K > 1, and we note that obtaining this result via combinatorial 
arguments seems challenging. Numerical approximations of the ODEs validate this asymptotic 
analysis. 

While the behavior of the threshold for the appearance of a giant component is similar for 
K 3> 1, combinatorial arguments yield that for all K > 0, t y (K) = 1, whereas t^(K) = max{^, K}, 
hence t^/t v = K for every K > I . 

It is possible to implement both processes Gr and Q y K efficiently by choosing an appropriate 
data-structure for holding the isolated and non- isolated vertices. Our implementation requires 0(n) 
memory and runs in time O(nlogn), and its results validate the above analytical results. 

1.2 Notations and main results 

A property of graphs is a collection of graphs closed under isomorphism. A property is said to 
be monotone (increasing) if it is closed under the addition of edges. Throughout the paper, we 
say that a property of graphs on n vertices occurs with high probability, or almost surely, or that 
almost every process Q satisfies this property, if the probability for the corresponding event tends 
to 1 as n — > co. Note that, when proving that certain statements hold with high probability, one 
may clearly condition on events which hold with high probability. 

On several occasions, we examine the processes Q y K or starting from some arbitrary graph 
H (instead of the edgeless graph). We denote these processes by Gk\h and Gk\h- 

Given a process G, we let G % denote G after % edges. It will be convenient to use two timescales 
when referring to G, which we denote by: 

g(t) -.= G tn/2 , 
g[t] -.= g*f lo s™. 

We say that t is a threshold for the appearance of a giant component in G if for every e > 0, 
with high probability G{t — e) does not contain a giant component and yet 0(t + e) does contain 
one. Similarly, we say that t is a threshold for connectivity if for every e > 0, with high probability 
Q[t — e] is disconnected whereas (/[i+e] is connected. According to these definitions, both thresholds 
equal 1 for the original random graph process Gi- 
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Theorem 11.11 determines the threshold for connectivity in both processes, t^ and t£ , and shows 
that = 1 for any K > whereas t£ = K for any i^T > i, hence their ratio is K for any X > |. 
Furthermore, as we later state, it follows from that a ratio of [~max{-^, X}] is the maximal 
possible between the threshold of and the threshold of Q\ for any monotone property, and 
therefore, £t achieves this maximum. 



Theorem 1.1. For every K > 0, t^(K) = 1 and yet t£(K) = max{^, K}. In the special case 

l 
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K = we have %(0) = t£(0) ~ 1 



For a given graph G = (V, E) on | V| = n vertices, we let C = C(G) denote the set of connected 
components of G, and for i £ N, we define Cj = Ci(G) to be the set of components of size i: 
Ci = {C £ C(G) : \C\ = i}. Whenever it is clear as to which graph process Q we are referring, we 
use the abbreviation C\ to denote Ci{Q l ). The fractions of vertices which belong to components of 
size 1 and 2 are defined as: 

1(G) = M , HG ) = 2 -M , 
n n 

and the susceptibility of G, S(G), is defined to be the average size of a connected component, 
averaged over all vertices: 

S(G) = -^2\C(v)\ = - £ |C7| 2 , 

vdV CeC(G) 

where C(v) denotes the connected component of v. The relation between S(G) and the existence 
of a giant component in G is immediate: if G contains a component of size an for some a > 0, 
then S(G) > a 2 n, and if S(G) > an then clearly there exists a component of size an. Therefore, 
G has a giant component iff S(G) = 0(n). 

The methods used in to analyze Q y K and describe t^ as a singularity point to a system of ODEs 
can in fact be applied to a wider class of graph processes. Namely, these methods can be applied to 
every process where the weight function W on the missing edges satisfies max Wf min W < K for 
some constant K > 0. Instead of repeating the complete set of arguments of in order to prove 
analogous results on Qj£, we summarize the arguments briefly, and proceed to prove the necessary 
conditions required for them to work. 

Following the ideas of ^ and previously of JO]) define the following system of coupled ODEs: 



y' 



-y 



l + (K-l)(l-y)2 j (1) 

2/(0) = 1 

z 2 + {K _ 1){z _ y) 2 

l + (K-l)(l-y)2 . (2) 

z(0) = 1 



We further define: 



w 



y 2 — 2wy — 2Kw(l — y) 



l + (K-l)(l-y)2 , ( 3 ) 

w(0) = 
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where y is the solution to (JTJ. Theorem 11.21 states that 1(G), S(G) and /2(G) along the process 
are approximated by the solutions to the ODEs above, and that tg is equal to the singularity 
point of the solution to (J2J) for any K > 0: 

Theorem 1.2. Let y(t), z(t) and w(t) denote the solutions for ([T ]) .p |) .((3* )) . and let x c denote the 
singularity point of z(t) if such exists, and 00 otherwise. For < 5 < 1, let t$ > be the minimal 
point satisfying y(r^) < 5. The following statements hold almost surely: 

1. For every 5 > 0, \I (G£(t)) - y(t)\ = o(l) and \I 2 (G&(t)) - w(t)\ = o(l) for all < t < r s . 

2. For every e, 5 > 0, \S (G^(t)) - z(t)\ = o(l) for all < t < mm{r s , x c - e}. 

3. For all K > 0, t£ = x c . 

The value of tg in the special case K = and the asymptotic behavior of tg are stated in the 
following theorem: 

Theorem 1.3. The threshold for the appearance of a giant component, tg(K), is a continuous 
function of K on (0, 00), and satisfies: 

where the o(\)-term tends to as K — > 00. 

The rest of this paper is organized as follows: in Section El we prove Theorem 11.11 by exam- 
ining second moments of processes which are easier to analyze, and provide results of computer 
simulations of t^ and t£. In Section |3J we prove Theorems 11.21 and II .31 using analysis of differential 
equations, and provide results of computer simulations of tg. 

2 Thresholds for connectivity 

The proof of Theorem 11.11 uses the following result of [Q : 

Definition. Let M S N. An M -bounded weighted graph process on n vertices, H = H(n), 
is an infinite sequence of graphs on n vertices, (Tt° ,T~C , ■ ■ ■), where TL° is some fixed initial graph, 
and Tl t is generated from H 1-1 by adding one edge at random, according to a distribution of the 
following type: the probability of adding the edge e to TL 1-1 is proportional to some weight function 
W t (e), satisfying: 

max WAe) < M min WAe) . 
If for some v > TL V = K n , we define TL 1 = TC U = K n for every t > v. 
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Theorem 2.1 Let H denote an M-bounded weighted graph process on n vertices, and let A 

denote a monotone increasing property of graphs on n vertices. The following statements hold for 
any t G N; 

Pr[W* G A] < Pr[Si M V G A , (4) 
Pr[£i V G A] < Pr[H Mt G A] . (5) 

Notice that for every > 0, the processes £/)^ and C?^ are both M-bounded, where M = 
[~max{ -j^ , K } . As it is well-known that I(t), the fraction of isolated vertices is (1 + o(l))exp(— C) 
at time Cn/2, Theorem 12 . II implies that the for any fixed C and K > 0, Q^{C) and Qj^{C) almost 
surely contain 0(n) isolated vertices. On the other hand, it is not difficult to see that Qq (1) has 
at most 1 isolated vertex, and we will later show that Qq (| + o(l)) almost surely has no isolated 
vertices. For this reason, we treat the cases K > and = separately. Lemmas 12.21 and 12.31 
below establish the precise behavior of these parameters when K > 0: 

Lemma 2.2. Let K > 0. For every < e < 1, £7#[1 + e] is almost surely connected, whereas 
G*k\\ — s] is almost surely disconnected. Altogether, = 1. 

Lemma 2.3. Let K > 0. For every < e < 1, £7#[(1 + e) max{^, X}] is almost surely connected, 
whereas Q^[K — e] and Gk\\ ~ £ ] are almost surely disconnected. Altogether, t£ = max{^,_RT}. 

As we stated in the introduction, combining the fact that t c = 1 for Q\ with Theorem 12.11 yields 
that t£(K) < |"max{p K}~\ , and indeed Lemma 12.31 shows that for K > 1, K G N, this maximum 
is achieved. 

The remaining case = is settled by the next lemma: 
Lemma 2.4. The threshold for connectivity for K = satisfies t^(0) = t£ (0) = ^. 



2.1 Proof of Lemma I2T21 

Let < s < 1; we first show that t^ < 1 + e. It is easy and well known that for every e' > 0, there 

exists some c = c(e') such that, with high probability, Qi{c) has a giant component of size at least 

el A 

(1 — e )n: assume that indeed this holds. Take < e < min{ - — , K} and c = c(e') as above, and 
notice that: 

1 + e> ~ 1 + e V ; 

By Theorem l2.ll the giant component of £/^ c n ^ 2 typically contains at least (1 — e')n vertices, where 
d = \m&x{K, ^}]c. Let H l denote the largest component of G^ 1 , and let A l v , v G V, denote the 
event: (v ^ H l ). The following holds for every i > dn/2: 

PrK4 i+1 l^l - ggl > K(1 " £ > > Izl (7) 



K(Q)-i)-(^-l)(l c 2 II)-^) + (l c 2 il)-(l + £ ')l 
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where the last inequality is by the fact that < e' < K. Setting To = c'n/2, this gives the 
following upper bound on for v ^ H T °: 



hence a union bound on the vertices of V \ H T ° implies that, with high probability, ^^-[1 + e] is 
connected. 

The lower bound on V c is slightly more delicate. A simple way to show the lower bound in Gi 
is to apply a second moment argument on the number of isolated vertices of However, in 

the case of Gk, using uniform upper and lower bounds on C\ (such as e'n and o(n)) does not yield 
useful bounds on the above second moment. The following claim resolves this difficulty : 

Claim 2.5. Let < 5 < \, and let {Hq} ti denote a family of arbitrary graphs on n vertices with 
(1 + o(l))ra 1-5 isolated vertices. If "H is a biased G\ process on n vertices which begins with Hq, 
that is: H ~ Gk\h , then TC[1 — 25] almost surely contains isolated vertices. 

Indeed, to obtain the lower bound on t y c from the above claim, fix 5 = e/2, and let r denote the 
minimal time at which < n 1-5 . Taking Hq = G^ T , Claim 031 implies that Gk\h almost surely 
has isolated vertices at time r + (1 — 25)^ logn, and in particular, we obtain that > 1 — e, as 
required. It remains to prove Claim 121)1 

Proof of claim. Let Hq and 5 be as above, and let Ti. ~ Gx\n a - Let B l u denote the event that the 
vertex u is isolated at time i, where i < nlogn. Clearly: 





(8) 



Thus, for T = (1 + e)^ log n we get: 



Pr[^] < exp 




PrhBj^lfl*] 



K{n-\e\\) + {\e\\-i) 




By the assumption on Hq, \C\\ < \Cl\ = (1 + o(l))n 1 s , thus: 




2 + 0{n- & ) 



n 



0) 



Similarly, we can define B l uv = B l u 



A B l v for i < nlogn and get: 




u,vl 



2K(n - \C\\) + 2{\C\\ - 2) + 1 _ 4 + 0(n~ s ) 



(10) 
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Let T = (1 — 26) T£ logn, and define Y = J2 U 1 b t = \Cf \ to be the number of isolated vertices of 
7i T . A straightforward second moment consideration implies that Y > almost surely. To see this, 
first note that @ along with the well known bound 1 — x > e~ x /Q~ x ) for < x < 1 yield: 

EY = |C?| p Prf^lK- 1 ] = (1 + o(l)y~ s (l - 2 + °^ 8) ) T > 



> (1 + o^n 1 - 6 exp I - v " ' ^gnF (l - 2<5)^ log n I = n i-* n -i+^+»W = n «+°W . (n 



(2 + o(l))/n _ n \ !_j 1+M+0(1) , +0(1) 



1 



2+o(l) 



By © and ifTUjl . there exist p = p(n) and q = q(n) such that Pr[-iS*+ 1 |S* J > p = 4+0 (" — - and 
Pr[-i£* +1 |.a£] < g = 2+Q ) t "~' ) . The following holds: 



Cov(1 b t,1 b t) = Pr[5^] -Pr[fl^ < (1 - (1 



< (1-25) log n-O^- 5 ) exp + o(l))-(l - 2<y)~logn) = n - 2 +M+o(i) _ 



Therefore, 

^ 2 Cov(1 b t, 1 b t) < n 2(l-«5 )rl -2 + 3 5+ o(l) = n «5 + o(l) = Q ((Ey)2) 

As Ey = w(l), applying Chebyshev's inequality gives: 

VarF Ey + V V Cov(1 b t , 1 b t) 1 

L J - (Ey) 2 ~ (Ey) 2 Ey w w 

This completes the proof of the claim and of Lemma 12.21 



2.2 Proof of Lemma EH 



Let < e < 1. The bounds t£ < (1 + e) max{±, if} and t£ > K - e will follow fr om arguments 
similar to the ones in the proof of Lemma 12.21 whereas the bound t£ > \ — e requires more work. 

To prove that t£ < (1 + e) max{^, K}, take s' > which satisfies: 

(12) 



[l-e') > 



1 + e 



1 



l-e / ) 2 + (l-^ / )^< + 



-A2 ; 



For instance, the following choice is suitable: 

' eK 



e < mm 



1 



1 — — minji^, 1} 
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As before, take c = c(e') such that, with high probability, Gi{c) has a giant component of size 
at least (1 — e')n. Defining H l to be the largest component of G^ 1 , Theorem 12.11 implies that 
\H c ' n l 2 \ > (1 — e')n with high probability, where d = [max{ if, j^}]c, and therefore we may assume 
that this indeed holds. Let A l u denote the event that a vertex u £ V \ H l , which belongs to some 
connected component C, joins the giant component at time i + The following holds: 

PrK] = mW °^ > where w c = { ' * ^ = 1 . (13) 

+ \K\C\ otherwise 

By the choice of c', for all i > c'n/2 we get: 

Prf.4M> w c {l-e')n 2w c {l-e') 



) + (if - 1) (-f i) n (i - (i - Jaip + if (i - M)3 



For if > 1, the denominator in the expression above is clearly bounded from above by nif , and for 
< K < 1 it is bounded from above by n (l — (1 — e') 2 + (1 — £ ; ) 2 if), hence (fH?)) implies: 

1 ui ~ nK{l + e/A) ~ nK{l + e) 1 ' 



Take T = c'n/2 and T = max{^, if }(l+e)^ log n, and define the following event for every u 

Bl = {ui H T " 



The definition of wc in f|13|) implies that > min{l,2if} and when combined with ()14j) this 
gives: 



^n K mm ™r /4> )^ 

/ l + e/4 mini i, 2} / 1 _ , « 
< exp — (max{-,if}(l + e)- logn - c - j j = n 



Therefore, the expected number of vertices of V\H T ° which do not join H T is clearly oil), implying 
that 1 + e) max{^, if }] is almost surely connected. 

The lower bound t£ > if — e follows from the following claim, analogous to Claim [231 

Claim 2.6. Let < 5 < |, and Zei {ifo}n denote a family of arbitrary graphs on n vertices with 
(1 + o(l))n 1_<5 isolated vertices. If Ti. ~ {?/^|# , W[if — 25] almost surely contains isolated 
vertices. 

Proof. As in the proof of Claim 12.51 let B l u denote the event {u is isolated at time i) and let 
B % uv =B l u A B*, where u,v £ V(H ). Clearly: 

ft \ 

Prh^Vi] ■ " 
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The assumption on Hq gives \C\\ < \C±\ = (1 + o(l))n 1 <5 , and thus for all i = O(relogn): 

[ " ' ^ ~ Kn (l + O(^) + O(n-)) " ' (1B) 

and: 

Pr h i^W ] = ^ = 4 + ° {n ~ S) (16) 

[ u ' J ©+(K-i)r^i)-^ [ ) 

Henceforth, a similar calculation to the one in the proof of Claim 12*31 gives the required result. ■ 



Choosing 5 = e/2 and applying the last claim implies that t£ > K — e, and it remains to show 
that t£ > g — s. This follows from the fact that Gx\h ~ £ ] contains components of size 2 almost 
surely. Unfortunately, we cannot repeat the last argument in order to show that indeed this is the 
case, since we have no guarantee that at time r, the hitting time for the property \C\\ < n 1_<5 , Gx 
still satisfies | C2 1 = re 1- " for some small a. To prove the next claim, which completes the proof of 
Lemma l2.2l we consider a simplified process, where it is possible to give a lower bound on IC2I, and 
show that this process stochastically dominates G^ with regards to |C2|. 

Claim 2.7. For all K £ (0, 1] and e > 0, G^i^ ~ £ ] almost surely contains a connected component 
of size 2. 

Proof of claim. Define the process G^, which is an approximated version of G^, as follows: at each 
step, assign a weight K to ordered pairs of the form {(u,v) : u, v ^ Ci}, and a weight 1 otherwise. 
If the ordered pair (u, v) chosen at some step i is a (self) loop or corresponds to an edge which 
already exists in Qki ^ ms s ^ e P is omitted. In other words, Gx assigns weights to all ordered pairs, 
and disregards selections of multiple edges or loops. Clearly, if G^ contains m < t edges, then 
G K ~ Gk" 1 - Furthermore, for any t = o(n 2 ) we have: 

E (* - \E(gg)\) < £ max{i K}^ = (o (±) + ° (~)) ' = °® > 

where E(H) denotes the set of edges of the graph H. In particular, for any fixed c, with high 
probability Gx( c ) ~ @k( c + (-Q) anc ^ @k[ c ] ~ ^k\ c + °(1)]> an d ^ * s sufficient to prove the claim 
for g*. 

Take < e < i. By well known results on the original Erdos-Renyi graph process, for any fixed 
t, with high probability: |C2(£7l(i))| = ©(re), with the constant tending to as t — > 00. For a short 
proof of this fact, one can verify that the following differential equation approximates the graph 
parameter I2 along Gi up to an o(l) error term: 

w'(t) = -2w(t) + y 2 (t) , w(0) = , 

where y(t) is the approximating function of the fraction of isolated vertices, I(t), and the time- 
scaling is of re/2 edges at each step. Substituting the well known fact (which also follows from 
Theorem 11.21 when substituting K = 1) that y(t) = exp(— t), it follows that w(t) = iexp(— 2t). 
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By Theorem 12.11 and the above fact, we deduce that for every sufficiently large c there exist 

< «i < a 2 < | such that Q^(c) almost surely satisfies a\n < IC2I < a 2 n. Let: 

eK 

c' = M{t > c : y(t) < — } , 

where y is the solution to the ODE Q), and set e' = y{d). By Theorem II, 21 with high probability 

1 [pK^>) = y( c ') + an< ^ m particular, Gj^(c') has (1 + o(l))e'n isolated vertices almost 
surely. Let Hq denote G^(c'), and take a\ = ai(c'), a 2 = «2(c') as above. According to these 
definitions, S = Ci(Hq) satisfies s = \S\ = (1 + o(l))e'n almost surely, and W = C 2 (Hq) satisfies 
a\n < \ W\ < a 2 n almost surely. Assume that indeed this holds. 

We consider the graph process Ti., which begins with Hq, and at each step selects an ordered 
pair (u, v) S V(Hq) 2 according to the following probabilities: 

■j^z if (u, v) is incident to W and S 

if (u, v) is incident to W and V \ S , 
A otherwise 

where the value of A > is chosen such that the probabilities sum up to 1. This is possible, since 
the probabilities for the first two types of pairs sum up to at most: 

2a 2 {e'/K + (1 - e')) < 3a 2 < 1 , 

as e' < y and a.2 < \- We claim that the following two events occur almost surely, and complete 
the proof of the claim: 

C 2 (H[l-s])nW^® , (17) 

\c 2 (H[\ - e]) nw\< \c 2 (g^\ Ho [^ - e]) n w\ . (18) 

A standard second moment consideration proves that (|17|) occurs with high probability. To see 
this, set: T = (^ — e) § logn, and notice that: 

2T 1 / e' \ 1 

(K(n - s - 2) + s) = -(1 - 2e) (1 - e') + — + o(l) logn = -(1 - e" + o(l)) log n (19) 



Kn 2 v v 7 7 2' 'V 7 2 

for some £ < e" < 2e (by our choice of e'). Next, let the random variable X e (e € W) be the 
indicator of the event: (e £ C 2 (H[^ — e])), and let X = ^ e£ w X e . ^he following holds: 

n*.=ii=(i 2A ' ( "~',;. 2)+2s ) T 

4if(n-s-2l -4ET + 4s x 2 



ifn 2 

Pr[X e = 1 A X e / = 1] = ( 1 



Kn 2 

for every e, e' £ W 7 ". Thus, a calculation similar to the one in the proof of Claim 1231 gives: 
EX > \W\ exp (-(2 - °(l))T K{n ~K~ 2 2)+S ) > \W\n(- 1+£ "V 2+ °W = , 
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and: 



Cov(X e ,X e ,) < 



Kn 2 



AK 



exp 




T < 4n- 3+e "+°( 1 ) . 



Therefore: 



Yl Cov(X e ,X e r) <4\W\ 



2 -3+e"+o(l) 



o(EX) 



e£W e'£W 



and in particular, Var(X) = (1 + o(l))EJ = o (EX) 2 ) and by Chebyshev's inequality we deduce 
that X > almost surely. 

It remains to prove that (|18|) occurs almost surely. This is achieved by a coupling argument: 
we claim that there exists a coupling of the processes G^\h and TC whose support consists of pairs 
(G t ,H t ) such that: G t ~ 0£|# , H t ~ U\ and (C 2 (fT t ) n W) C (C 2 (G t ) n W)- This clearly holds 
for t = 0, and by induction, it remains to extend the coupling from (Gt,H t ) to (C?t+i, -Ht+i). For 
this purpose, we apply the following lemma of (Lemma 2.2), which was first proved by Strassen 
in a slightly different setting: 

Lemma 2.8 Let U, V be two finite sets, and let R C U xV denote a relation on U, V . Let \x 

and v denote probability measures on U and V respectively, such that the following inequality holds 
for every A C U: 



Then there exists a coupling ip of [i,v whose support is contained in R. 

Let U, V be the set of all n 2 ordered pairs selected, representing the next pair selected by Gt 
and Ht respectively. Let [i denote the probability measure of each selection in U by Gt, and let v 
denote the probability of each selection in V by Ht- Define: 



By the induction hypothesis, Y C. X, and we define R to be {((u,v), (x,y)) : Y XtV C X u>v }. 

Clearly, if (u,v) is a loop or an edge which already belongs to Gt, it has no effect on X, and 
as Y x y C Y C X we get (u, v)R(x, y) for all x, y. Furthermore, every (u, v) which is not incident 
to any e £ Y also satisfies (u, v)R(x, y) for all x, y, as the components (u, v) may remove from X 
already do not belong to Y. Therefore, (|2Uj) is satisfied for every A C U which contains such pairs 
(u,v). 

It remains to prove that (|2T)|) holds for sets A C U such that A n E(Gt) = 0, and A consists 
entirely of edges incident to edges of Y. Let A = {e±, . . . ,e m }, and notice that ej ^ E(H t ) for 
all i, as Si is incident to some component e G Y of size 2 in Ht, satisfying e ^ ej as 7 C E(Gt)- 
Furthermore, for all i we have eiRei, since all components that removes from X in X ei are also 
removed from Y in Y ei . Thus, showing that yu(ej) < z^(ej) for all i will imply that A satisfies the 
condition of (|20j). Indeed, if = (u,f) is such that u,v ^ S (u,v are both non-isolated in Hq), 



fJ,(A) < v{{y £ V : xRy for some x £ A}) . 



(20) 



x = c 2 (G t )nw 

Y = C 2 (H t )nW 



x u , v = c 2 (G t u («,«)) nw 
y x , y = c 2 (F t u(x, y ))niy 
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then: 

K K 
M(ei} " n2 + (^-l)(n-d(G*)) 2 ~Krt~ ^ ' 

and otherwise: 

1 1 

K€l) ~ n2 + (^-l)(n-d(G*)) 2 " tfn* " ^^i) , 
by definition of the process TC, completing the proof of the claim and the proof of Lemma 12.21 



2.3 Proof of Lemma 12.41 

To prove that (0) = \ } we recall the following easy facts stated in pQ: by definition, Qq adds 
edges between pairs of isolated vertices until no such pair is left. Hence, after adding \n/2\ edges 
there is at most 1 isolated vertex in Q^, and Q y K behaves as Q\ on [n/2\ components of size 2 (and 
possibly 1 additional isolated vertex). Thus, tg (0) = 1 + % = §• 

Proving a connectivity threshold of + o(l))^ logn, we may assume n is even: otherwise, the 
single isolated vertex from time \n/2\ becomes connected almost surely at time co(n), and its edge 
set accounts to at most n — 1 = o(n log n) edges, not affecting the threshold for connectivity. By 
the discussion above, for even values of n, Qq has a linear number of components of size 2 and 
no isolated vertices at the time of appearance of the giant component. Thus, we deduce from the 
arguments used for the proofs of Lemmas 12.21 and 12.31 that, with high probability, Qq[\ — s] still 
contains components of size 2, whereas Qq[\ + e] is connected. Altogether, (0) = i. 

It remains to show that t£ (0) = \. Substituting K = in equation (JJ), it reduces to the form: 

✓ = ^,y(0) = l, 

provided that y ^ 0. This yields the solution: 

y{t) = 2- VT+2i . (21) 

Notice that y(t) is strictly monotone decreasing, and reaches at t = |. Therefore, if we denote by 
tq(Qq) the minimal time t at which I(Go(t)) = 0, Theorem 11.21 implies that for every < S < 1, 
|r — || < S almost surely. To see this, let xq be such that y(xo) = 5, and let x\ = max{xo, | — S}. 
By Theorem 11.21 with high probability: 

— S))>I (g A (x!)) = y{x x ) + o(l) > \y{x x ) > , 

where the second from last inequality holds for sufficiently large values of n. On the other hand, 
y{x\) < 5, and hence I (Qq (| + 5)) < (£?o ( x i + <^)) = almost surely, since each edge eliminates 
at least one isolated vertex in Qq . 

Equation © takes the following form after substituting K = and the value of y(t) as it is 
given in (|2"T|) . provided that y / 0: 

, , . y — 2w 1 — w 

w '(t) = * = 2 - - 1 , w(0) = , 

W 2-1/ ^/TT2t 7 K J 
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Figure 1: Comparison of the numerical results for and t^, and estimations of t£ and according 
to computer simulations of the model. 



which has the solution: 



w(t) 



1 



3 2(l-^l+2t) _ i^2t . 



(22) 



4 4 2 

~ 0.1485, hence from the discussion above we obtain that 



For t = | we have w(t) - ( 
# = Gk(t) almost surely satisfies h(H) = Q(n) (and 1(H) = 0). From that point on, the process 
is equivalent to Gi\h, and from the similar arguments to those used in the proofs of Lemmas 
l2~2l and E3] we obtain that % (0) = \ . ■ 



2.4 Computer experiments of t£ and t^ 

Maintaining the set of isolated vertices and the edges already added to the process provides all the 
information needed to add the next edge to and Q y K at an O(l) cost. In order to recognize the 
threshold for connectivity, the set of connected components must be efficiently maintained. Our 
implementation holds the components in linked-lists, according to the Weighted-Union Heuristic 
(see, e.g., jSj p. 445) which guarantees an average cost of O(logn) for uniting components. 

Figure ^ shows the results of t£ and according to simulations of both models on n = 10 4 
vertices. The values of t£ and were averaged over 100 tests per value of K. 
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3 The appearance of a giant component in 
3.1 Proof of Theorem fTT2l 



We begin with a short summary of the methods used in to analyze t^(K) for K > 0. First, the 
authors prove Theorem 12 , 1 1 and deduce that the process Q y K is stochastically dominated by Q\, up 
to a timescale factor of |"max{^, K}~\ . By applying the differential equation method of Wormald 
|12| to the approximated process Q y K (which selects an ordered pair at each step, similar to the 
process introduced in the proof of Claim l2~Tj) . the parameters 1(G) and S(G) are approximated 
by y(t) and z(t), solutions to a system of coupled ODEs. Henceforth, a repeated use of Theorem 
3.1 of ^Uj, which relates the susceptibility and the appearance of a giant component, implies that 
the singularity point z(t) is equal to t y c . 

We note that the methods of can be applied to any family of M-bounded processes Vm ( as 
referred to in Theorem 12.1(1 provided that the following conditions hold: 

1 . Attempting to approximate a fixed number of bounded graph parameters Y\ , . . . , (such as 
1(G) and h(G)) by functions yi, . . . , yd for 1 < t < T, we require that 

E «VW = I (y l( ^ )r .. )y ^ ))t ) + erryi(t) , (23 ) 

for all i and t, where max^j err yi (t) = o(l). 

2. In addition to the above approximations of Y\,...,Yd by yi,...,yd, when attempting to 
approximate S(G) by the function z(t) for 1 < t < T, we require that 

E ^P^_S(n^ = | (Y^l.-.^^S^t) +err z (t) , (24) 
where maxterr^i) = o(l). 

To show the above, we may use the fact that with high probability, the largest component is 
of size O(logn) as long as t < x c — e, where x c is a possible singularity point of z(t). This 
follows from the proof of Theorem 1.3 of [Q, which applies to this generalized setting as well. 

3. Finally, if the above function z(t) has a singularity point, we require that it is uniformly 
bounded (regardless of M). 

If the above 3 conditions hold, we obtain that Y\, . . . , Y^, S are within o(l) distance from y\, . . . , y^, z 
along the process Vm for 1 < t < T. Furthermore, if z has a singularity point at x c and the above 
conditions hold with T = x c — e for any e > 0, it follows that the appearance of a giant component 
in Vm is at i = x c . 

We first show that the above 3 conditions hold for K > under the assumptions of Theorem 
11.21 By the well known fact that for any constant t, I(Q\(t)) = Q(n) almost surely, Theorem 12.11 
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implies that ts can be taken to be arbitrarily large. In particular, we can take t§ > x c . In order 
to verify that (|23|) holds for /(G) and 12(G), let G and G' denote Q 1 ^ 1 and Qj^ t+l respectively for 
some fixed t, and let I = 1(G) and I2 = 12(G). We have: 



and: 



E 



2/n v ; l + (if-l)(l-/) 2 v ; l + (if -l)(l-i) 2 

= l + (K-"lKl-/) 2 +0 ^' <*> 

h{G')-h{G) _ 1(1 - |) - 2J 2 J - 2ifi 2 (l -I-h)- 2KI 2 (h - f) _ 
2/n ~ 1 + (if — 1)(1 — I) 2 

- ^r^z^^ ^W")^^)). (26) 

where in both cases we used the fact that if > to obtain an upper bound of o(l) on the error 
term. To prove (|24j). set 5 = S*(G) and observe that: 

, 5(0') -S(G) 



E- 



2/n 



™/2 ^i(i-i)- + 2i y M 2 M+ y y tficilic^idiic-ai 

T7. -n ^— T7. n -J -J 



1 + (if — — I) 2 \ n n ' n n ' n 2 n 

|C|>1 |Ci|>l |C 2 |>1 

27(5 - /) + i 2 + if (S - I) 2 



+ 0(//n) +0 (£^) 



l + (^-D(l-/) 2 

|C|>1 



5 2 + (if- i)(5-i) 2 \A.n( l C H ^ 

|C|>1 



l + (if 



and the assumption that \C\ = O(logn) gives a bound of o(l) on err z . The next claim therefore 
completes the statements of Theorem 11.21 for the case if > 0: 

Claim 3.1. For every if > 0, tg(K) < 5. 

Proof. Consider the case if > 1, and let y(t) and z(i) denote the solutions to the ODEs © and 
<PJ) respectively. Recalling that y(t) < 1 for every t > 0. <P)l yields that z'(t) > provided that 
<z(i) > 1- Thus, the initial condition z(0) = 1 implies that z(t) is monotone increasing in t, and in 
particular: 

, _ z 2 + (if - l)(z - y) 2 ^ z 2 + (K- l)z 2 (l - y) 2 _ 2 

l + (if-l)(l-y) 2 " l + (if-l)(l-y) 2 ' {6) 

where inequality is by the fact that y > and z > 1. By standard considerations in differential 
analysis , (|28|) and the initial condition z(0) = 1 imply that z(t) > for every t > (as = j—l 
satisfies z! = z 2 and z(0) = 1), and in particular, tg < 1 for any if > 1. 
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We are left with the case < K < 1. Clearly, (JJ) implies that y' < for every t > 0, and 
furthermore, y' < as long as y > 0, hence y is strictly monotone decreasing from 1 to 0. Let t* 
be such that 



y(t*) = 1 - VT^K . (29) 



Notice that the solutions to the equation —x 2 + 2x — j^f^ = are: x\p = 1 ± y 1 — if 
< K < |, and no solution exists if | < < 1. In both cases, — x 2 + 2x — < for every 
x < 1 — \/l — -K"- Therefore, the definition of i* and the fact that y is monotone decreasing give: 

-y(t) 2 + 2y(t) - — - < for every t > t* , 

or equivalently: 

(1 - K)y{t){2 - y{t)) < K for every t > t* . (30) 

Rewriting equation © as: 

/ _ z 2 _ (1 _ g)^ _ y )2 _ ^2 + (1 _ ^ )y(2z _ y) 
Z l_(l_^)(l_y)2 ^ +(1 _^ )2/(2 _ y) • I ^ 

it follows that for every t > t* , z' > \z 2 . Furthermore, for every t > 0, z 1 > 0, and hence 
z(t*) > z(0) = 1. We obtain that the function w(t) = z(t — t*) satisfies w' > \w 2 for every t > 
and w(0) > 1, and by the same consideration as above, w(t) > for every t > 0. Altogether, we 
deduce that < £* + 2, and it remains to provide an upper bound on t* . 



2' 

■t lilll //'( / \ — - I I - — / '■" 1 C M I K' i U K »/;' "> _ ( 

2-t 

< / ; 4- 9 Mnrl il i'ciiV.iiiK tri nrrivnrlo mii mvrw>r lin 

For this purpose, define u(t) = 1 — y(t), and consider Q for < i < i*: 

. . 1 — u 1 — it 1 

M = -V = -i T, 77\ 9 > 



1 - (1 - K)u 2 - 1 - it 4 1 + it + it 2 + it 3 ' 

where the inequality is by the fact that u{t) < y/1 — K for < t < t* . Define w(t) to be the 
solution to the differential equation: 

w ' = n / 2 1 3 > ™(°) = . ( 32 ) 

it follows from the above mentioned argument that u(t) > io(i) for < £ < £* . The solution to (|32|) 
satisfies: t = Y^j=i hence w(to) = y/l — K for to = Ylj=i ~ ~p — — f§- As M (*o) > ^(^o)> it 
follows that t* < to < || , completing the proof. ■ 

In the special case K = 0, Qq is no longer an M-bounded process, however the assertions of 
statements 1, 2 of Theorem II . 21 remain valid and follow from (|23f) and (|24|). by applying Wormald's 
differential equation method directly. To see that Q23[) holds, notice that as long as I(C/^-(t)) > 5 
for some fixed <5 > 0, the denominators in (|25j) and (|26|) are 0(1), and the approximation remains 
valid (note that z(t) has no singularity point for K = 0). To show that (j2H> holds, set S = S(G), 
and note that: 
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Figure 2: Comparison of the numerical results for tg, and estimations of tg according to computer 
simulations of the model. 



E 



S(G') - S{G) 
2/n 



n/2 



[l -I) 



2 » ran ^ 

C6C 

|C|>1 



n n 



2I(S-I)+I 2 _, r/ . 25-/ _, T . , 
+ 0(J/n) = - — - + 0(J/n) 



21 -P 



I 



(33) 



where we used the fact that I > for t < r^. This implies that the error term err z is o(l) without 
making any assumptions on the size of the largest component, and completes the proof of the 
theorem. ■ 



3.2 Computer experiments of tg 

We conducted simulations of tg using the implementation of mentioned in 12.41 In these sim- 
ulations, the number of vertices was n = 10 6 , and the threshold for the appearance of the giant 
component was taken to be the minimal time at which contains a component of size an, where 
a = 0.01. The value of tg (K) was averaged over 10 tests for each value of K. 

Figure |2] shows the comparison between the values of tg according to the above computer 
simulations, and the values obtained by numerically solving the ODEs Q and (J2J by Mathematica. 
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3.3 Proof of Theorem fOl 



First, the fact that t^(K) is continuous follows from the general continuous dependence of ODEs 
on their parameters (in this case, the single parameter K). 

For the special case K = 0, recall that in our treatment of Gq for the proof of Lemma 12.41 we 
showed that for every 5 > 0, with high probability y(t) = 2 — ^1 + 2t approximates I(Gq (t)) for 
t < | — S, and I(Gq (§ + 5)) = 0. Take 5 > and consider the interval [0, | — 5]; equation (J2J) takes 
the following form when substituting K = and the solution to y{t): 

= 2?-(z- yf = 2z-y = 2{z - 1) 

22/-y 2 2-y ^/T+2~t ' V j 

Taking w = z — 1, we obtain the linear equation: 

u/ ; 2 w = 1 . (35) 

VT+2i 

Multiplying 1)35|) by its integrating factor and integrating by parts, we get: 



w(t) = exp |2 j(l + 2t)~^dt^ J exp j-2 ^ (1 + 2t)~^dt 



dx 



~{VT+2i+ i) + Cexp(2VTT2t) . 



The initial condition u)(0) = gives C = |e 2 , hence: 



Z ( t ) = L KVT+Tt-l) _ I^T+2t + ? for all t € [0, - - S] . 
4 2 4 2 

The above solutions for y(t) and ,z(t) give y(|) = and 

*(3/2) = ^e 2 - i . (36) 

Let To = to(^q ) denote the first time t at which <5q * has no isolated vertices, i.e., the hitting time 
for the property: {G : 1(G) = 0}. By the above arguments, we have: 

r = (l + o(l))n/2 (37) 



2 

almost surely. We claim that proving that with high probability: 

S(C?n = z(3/2)+ (l) (38) 

implies the required result on tg (0). Indeed, once no isolated vertices are left, the process Go adds 
edges according to the uniform distribution, and hence from that point the susceptibility follows 
the equation z' = z 2 (e.g., see the case K = 1 of the analysis of Gx or Gjc)- By (|37|) . we obtain 
that for some x c > |: 

S(G£(t)) = (1 + o(l))^— t for any | < f < x c , 
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and the value of x c is derived from the initial condition ()36l) : 

W= x c = l + ^i- ( 39 ) 

It is left to show that (|38jl indeed holds. The lower bound S(Qq T °) > -z(l) — o(l) follows from 
Theorem 11,21 which states that z approximates S until time | — S for any (5 > 0, and from the 
continuity and monotonicity of z. That is, for any fixed £ > 0, choosing a sufficiently small (5 > 
such that z(| — (5) > z(|) — £ gives: 

5 (s A (^ - *)) = z{\ -5) + > Z (|) - £ + o(l) . 

For the upper bound S(£/q T °) < — o(l) we are required to examine the second moment of 
S (Qq (| + <5))- Assume by contradiction that: 



Pr 



■2 

and choose 5 > small enough such that: 



> a for some fixed a, £ > , (40) 



< a , and: (e 4<5 - l) Z (|) < £/2 . (41) 



Set To = (| — <5) | and T\ = tq — 1, and note that, with high probability, A := r ^J° satisfies 
5/2 < A < 25, and therefore we may assume that this holds. We consider Qq for T = Tq, . . . ,T%, 
and let S = S (Gq T °) ■ As I(Gq T ) > 0, the calculation which yielded (jS) gives: 



and hence: 



Therefore: 



i + ^ s(s aT ) < e (s(s AT+1 ) | s(g£ 

e A S < ES(G£ T0 ) < e 2A S . 



Similarly, we can write the expression for the second moment of S(Qq 



e(^o aT+1 ) 2 I s(a AT ) 2 ) <s(f# T ) a [i + n o _ jiiyT] 



2-/(S A 
and hence: 

E(S(£ AT0 ) 2 ) <e 6A 5 2 . 

We obtain that: 

VarS(£ AT °) < (e 6A - e A ) 5 2 < (e 125 - e 2S ) S 2 , 



20 



500 1000 



5000 10000 



_, 1 — k 

50000 100000. 



Figure 3: Comparison of the values of tg (K), obtained by numerical solutions of the ODEs 
and the asymptotic approximation of Theorem II .31 Logarithmic scale was used in both axes. 



where the last inequality is by the fact that A < 25 almost surely. By ()41|) and the fact that 
So < z(|) almost surely, this implies that: 



\^S{G^)<a(i/2y 



(42) 



By (|41j) we have: 



Pr 



< Pr 



< 



<Pt[\S{G, 



Ann 
J 



m 



A TO I 



> 



and thus combining Chebyshev's inequality with (|42|) gives: 



Pr 



< a 



contradicting the assumption (|40|). 

It remains to show that for K S> 1, tg(K) = (1 + °(1))^7| fl + fj) Consider equation 
the initial condition y(0) = 1 suggests that we examine the function u(t) = 1 — y(t), which 
satisfies u(t) <C 1 near the origin: 

, u — 1 



-u 



u 2 {K-l) + l ' 
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and substituting the initial condition of y we have: 



As u <C 1, the w 2 -term at the denominator is negligible, and we have: 

u'^^ TT ,u(0) = 0, (44) 

where here and in what the follows, we denote by ~ equality up to leading order terms. Rearranging 
the last equation to the form Ku 2 u' + v! ~ 1 and integrating, we obtain that, up to leading order 
terms, u satisfies the equation: 

— u 6 + uf*t. (45) 

We note that this immediately gives u ~ t for u -C H=, however we are interested in the behavior 
of u precisely when t = 0(— L). Applying Caradano's solution to the above cubic equation gives 
the following approximation of u when t <C 1 (and hence u <C 1): 

1/3 / , \ 1/3 



3t / 1 / 3t \ 2 \ / 3t / 1 / 3t x 2 



Moving on to z(t), we substitute u> = z — 1 in equation @ and obtain the following: 
, (w + l) 2 + (K - l)(w + u) 2 (w + l) 2 + (K - l)(w + u) 2 



1 + ifu 2 - u 2 1 + Kn 2 



Next, we may replace w + 1 by 1 whenever toil, and furthermore, whenever u; = £1(1) clearly 
the dominant term is K(w + u) 2 . Altogether, we obtain the following uniform approximation for 
w: 

K(w + u) 2 + l , 2 , , 

W Ku 2 + 1 ^ {W + U) + > U ' 

Adding u' to both sides of the equation and rearranging, we obtain: 

w' + v! 
K(w + u) 2 + 2 

hence if we define v = \J~^(w + u) we obtain: 

v* + 1 

and thus: 

v ~ tan(V2Xn) . 

Returning to z, we get: 

z«l| J~^tan(V2Ku) — u . 
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This implies that the singularity point x c satisfies y/2Ku(x c ) = ^. Recalling equation (|45l) . we 
have: 

K . .o . K 7T 3 7T 7T / 7T 2 \1 

x c — n(x c ) + u(x c ) = 7 — H - = — = 1 + 



3 v ' K ' 3 16KV2K 2V2K 2^2 \ 24 J ' 
Figure El shows the excellent agreement between the above asymptotic approximation of tg(K), 
and its value as obtained by numerically solving the ODEs (^Q) and (j2J by Mathematica. ■ 
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